Cotunneling at resonance for the single-electron transistor 
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We study electron transport through a small metallic is- 
land in the perturbative regime. Using a recently developed 
diagrammatic technique, we calculate the occupation of the 
island as well as the conductance through the transistor in 
forth order in the tunneling matrix elements, a process re- 
ferred to as cotunneling. Our formulation does not require 
the introduction of a cut-off. At resonance we find significant 
modifications of previous theories and good agreement with 
recent experiments. 

Electron transport through small metallic islands is 
strongly influenced by the charging energy associated 
with low capacitance of the junctions pHq]. A variety of 
single-electron effects, including Coulomb blockade phe- 
nomena and gate-voltage dependent oscillations of the 
conductance, have been observed. Usually, they are de- 
scribed within the "orthodox theory" |jj which treats 
tunneling in lowest order perturbation theory (golden 
rule) and corresponds to the classical picture of inco- 
herent tunneling processes (sequential tunneling). As a 
necessary condition one needs weak tunneling, i.e., the 
conductance of the barriers has to be low 



a = h/{4ir 2 e 2 R T ) < 1. 



(1) 



Despite the success of this straightforward approach, it 
was found experimentally and theoretically that there are 
several regimes where higher order tunneling processes 
have to be taken into account. 

First, in the Coulomb blockade regime, sequential tun- 
neling is exponentially suppressed. The leading contri- 
bution to the current is a second-order process in which 
electrons tunnel via a virtual state of the island. Averin 
and Nazarov M evaluated the transition rate of this "in- 
elastic cotunneling" process at zero temperature. At fi- 
nite temperature, divergences arise, but the authors of 
Ref. B provide an approximation which is valid far away 
from the resonances. They assumed that some regular- 
ization procedure will overcome the divergences. Their 
results were confirmed experimentally pj . 

Second, it was found recently ||f7| that even at res- 
onance, where sequential tunneling is not suppressed, 
higher order processes are important and can lead to 
a significant change of the conductance. Similar ef- 
fects were discussed for the average charge of the single- 
electron box in equilibrium ||- |Dl|l5[ . A diagrammatic 
real-time technique was developed for metallic islands 
||[7| as well as for quantum dots fL2| , [l3[ in order to give a 
systematic description of the various tunneling processes. 



The effects from quantum fluctuations were shown to be- 
come observable either for strong tunneling ao ~ 1 or 
at low enough temperatures aohiEc/T ~ 1, where Eq 
denotes the charging energy. The predicted broadening 
of the conductance peak as well as the reduction of its 
height was confirmed in the experiments of Joyez et al. 
Jl4| in the strong tunneling regime. Within the theory, 
only processes where the two classically occupied charge 
states are involved (even virtually) were included. There- 
fore, it was necessary to introduce a band-width cut-off 
~ Eq, which prohibits a comparison with experiment 
without fitting parameters. However, at low tempera- 
tures, a quantitative fit between theory and experiment 
is possible provided a renormalized value for the charging 
energy is used as it was determined in the experiment. 

In this Letter, we use the same diagrammatic tech- 
nique to obtain the total current in second order in ao 
including all relevant processes such that no cut-off re- 
mains. All terms are regularized in a natural way. This is 
important for a comparison with experiments since only 
system parameters enter the result, while a cut-off or 
some arbitrary regularization scheme would be an un- 
determined ingredient in the theory. At resonance we 
obtain new contributions compared to the existing the- 
ory of electron cotunneling. They emerge from a change 
of the occupation probabilities and a renormalization of 
the charge excitation energy. For realistic parameters 



T/E c 



0.05 and a. 



L — „,R 



0.02 the corrections are of 



order 20%. We compare with recent experiments |L4| and 
find reasonable agreement without fitting any parameter. 
The single-electron transistor is modeled by the stan- 
dard tunneling Hamiltonian 



H = H L + H I< +H I + H I 



ch 



Ht — Hn 



Hi 



(2) 



Here H x = £ A 



kn u ikn 



a rkn and Hj = J2 



qn e qnC qn C qn 

describe the noninteracting electrons in the two leads 
r=L,R and on the island where n is the transverse chan- 
nel index which includes the spin. The wave vectors 
k and q numerate the states of the electrons for fixed 
r and n. In the following, we consider "wide" metal- 
lic junctions with N ^> 1 transverse channels. The 
Coulomb interaction of the electrons on the island is 
modeled by H^ = Ec{n — « x ) 2 , where Eq — §q is 
the charging energy, and en x = ClVl + CrVr + C g V g , 
with C — Cl + Cr + C g , describes an external charge, 
which depends on the voltages and capacitances of the 
circuit. The excess particle number operator on the is- 
land is given by h. The charge transfer processes due to 
tunneling are described by 



tfn 



r=L,R kqn 



\*kq a rkn C( l ne 



-%(p 



h.c. 
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where T£™ are the tunneling matrix elements and e v 
changes the excess particle number on the island by ±1. 
The tunneling matrix elements T£™ = T r ™ are considered 
independent of the states k and q. They are related to the 
tunneling resistances Rt,t of the left and right junction 
via jji^ = ^- J2 n A r r(0)7V; l (0)|T rn | 2 , where Aj Tl (0) and 
-/V r "(0) are the density of states of the island and the leads. 
In the following we use the diagrammatic technique de- 
veloped in Ref. |6|,[7| . The noncquilibrium time evolution 
of the charge degrees of freedom on the island is described 
by its reduced density matrix, which we obtain in an ex- 
pansion in Ht. The reservoirs are assumed to remain in 
thermal equilibrium (with electro-chemical potential /i r ) 
and are traced out by using Wick's theorem, such that 
the Fermion operators are contracted in pairs. For a large 
number of transverse channels N, only those configura- 
tions contribute where the operators a rkn Cq n at one time, 
i.e., from one term Ht are contracted with the operators 
c qn a rkn from one other term Ht, and not from different 
ones. These "simple loops" dominate over more compli- 
cated configurations with more than two times connected 
by contractions. Matrix elements of the reduced density 
operator are visualized in Fig. (fil). The forward and the 
backward propagator (Keldysh contour) are coupled by 
"tunneling lines" (simple loops) associated with the junc- 
tions to the reservoirs r. Each tunneling line with energy 
u> represents the rate af 1 (oo) if the line is directed back- 
ward (forward) with respect to the closed time path with 



(3) tunneling) it reads p„a + (A„) — p) n ^_ 1 a^(A n ) — 0. At 
low temperature at most two charge states (n — 0,1) are 
important, all other states are suppressed exponentially. 

Due to higher order processes, however, the occupation 
is modified and also the probability for the other charge 
states can be nonzero (they are algebraically suppressed, 
but not exponentially). The Master equation expanded 
in second order gives a relation between the rates in sec- 
ond order £«,„' (diagrams with two lines) and the occu- 
pation in first order p n which lead to a deviation in the 
average occupation (n) = ^ n n Pn = ( n ) + (n)^ + . . .. 

The stationary current I r = -ie Yl n ,n' Pn^^ n ' 
through reservoir r uses the rates £„ „/ where the right- 
most tunneling line corresponds to reservoir r and is an 
outgoing (incoming) one if the rightmost vertex lies on 
the upper (lower) propagator (and vice versa for Ti r ~ n , ) . 
There are two types of diagrams contributing to the sec- 
ond order correction of the current 1^ : those of the form 
p(o) vj(2) anc j ^g others like p^E' 1 ' . The first ones corre- 
spond to "new" , second order processes, and the second 
ones are responsible for a modification of the first order 
processes due to the fact, that the occupation probabil- 
ities are changed in higher orders. The latter have not 
been considered in previous theories. As we will see later, 
both contributions are important. 

In lowest order the average occupation 



: M-±«o 



exp[±/3(w - /Mr)] - 1 



(4) 



They are associated with changes of the charge state, as 
indicated on the closed time path. Finally, we associate 
with each tunneling vertex at time t a factor exp (iAEt) 
where AE is the difference of out- and incoming energies. 
If the vertex lies on the backward propagator it acquires a 
factor —1. We define ao = JZ r ajj an( i oc{u) = ^ r a r (oo). 
The transport properties are affected by the energy 
difference of adjacent charge states 

A n = E ch (n + 1) - E ch (n) =E C [1 + 2(n - n x )] . (5) 

The formally_exact Master equation reads in the sta- 
tionary case 



Pn 



/ [Pn' 2-Hi' ,n Priori,'. 

n'^n 







(6) 



with the probability p n to be in charge state n and transi- 
tion rates £ n ,n' to change from state n to n' . In the per- 
turbative regime we write p n = Pn +Pn +Pn + ■ ■ ■ and 
=n,»' = ^1 + ^+^1 + . . . where p { n k) and E<£, de- 
notes the term ~ a§ of the expansion. The Master equa- 
tion must hold in each order. In lowest order (sequential 



\(0) 



a +(A )/a(A ) 



(7) 



is only smeared by temperature. Quantum fluctuations, 
however, yield 



\(i) 



_L d_ 

ZEcdn* 



^ 0) (Re f T-(A_ 1 )-Re ( 7 + (A )) 
+^ 0) (Re ( r-(A )-Rea+(A 1 ))l (8) 



where Rea ± (uj) = Er^LR^^M an d 



oo 

Reat(u>) = Re / duo 



uo-uo' + i0+ 



a 1 ± j U - ( W - Mr) 



ln^--Re* (i-£-(u>-Hr) 



(9) 



The integrals a (oo) are divergent since the integrand 
does not decay to zero for |w| — * oo. Therefore, we intro- 
duce a Lorentzian cut-off with cut-off parameter U. For 
the physical quantities like (n), however, only combina- 
tion of terms occur where this cut-off drops out, i.e., the 
divergences of different integrals cancel. 

In equilibrium, i.e. at V — 0, the transistor is equiva- 
lent to the single electron box. A systematic perturbative 
expansion of the partition function (up to order a§) was 
performed by Grabert jlq ]. The result Eq. rt|) is identi- 
cal to Grabert 's result in order ao, which, at zero tem- 
perature, reads (n)^ = aoln[(l + 2n x )/(l — 2n x )]. As a 



generalization Eq. QSj) also applies for the nonequilibrium 
situation, i.e. V ^ 0. 

The current I = Jl = — Ir is in lowest order given by 



/W(Ao) 



47r 2 ea L (Ao)a R (A ) 
h a(Ao) 



[/ L (A )-/ R (A )]. (10) 



The cotunneling contribution can be divided into three 
parts /( 2 )(A ) = YLi 7 f } (A ) with 



4 2) (A ) = - f duo iW (w)a(w) 



Re 



[p<°> ( x 


1 A 


P[) \uj-A + M+ 


w-A_i +i0+j 

1 VI 


" x v^- A i + *o+ 


w- Ao + iO+y 



(11) 



dujRe 



a (a;) 



/f ) (Ao)=I«(Ao) 



a (w) 



w-A o -M0+ cj-A_i+iO+ 



a+(w) 



a+(oj) 



uj - Ai + i0+ u - A + iO H 



(12) 



/ 



dwRe 



4 2) (Ao) = 



/ (1) (Ao) 



_9_ 

9A 

a _ (w) 



oj-A +iO+ w-A_i+iO+ 

a + (a;) a + (aj) 

w - Ai + i0+ uj - A + i0+ 



(13) 



The poles at uj = A are regularized in a natural way 
(it comes out of our theory and is not added by hand) 
as Cauchy's principal values Re , i0 + = P- and their 



derivative Re 7 



dx x ' 



S°) 



' (x+M+) 2 

Deep in the Coulomb blockade regime, we have p " 
1, Pi =0 and I^(Aq) = 0. Consequently, the first 

(2) 

line of I{ is the only contributing one. At T = 0, the 
integrand is zero at the poles, and we can omit the term 
+i0 + . This gives the well-known result of inelastic co- 
tunneling H. At finite temperature, however, the reg- 
ularization scheme is needed which is not provided by 
previous theories within second order perturbation the- 
ory Jig . Our result is also well-defined for T/0. 

Furthermore, we are able to describe the system at res- 

(2) (2) 

onance. In this regime, /j an( l ^3 become important. 
The origin of the second term may intuitively be inter- 
preted as the reduction of the first order contribution 
/^'(Aq) since quantum fluctuations lead to an occupa- 
tion of the adjacent charge states n = — 1 and 2 which 
is only algebraically suppressed, and no more exponen- 
tially. Therefore, the probability of the system to be in 



state n = or 1 (which is necessary for the first order 
process) is decreased. The third term may indicate the 
appearance of a renormalization of the excitation energy 
Ao J6p]9|,|ll[ . Due to this renormalization (which should 
be proportional to the gap and reduce it) the system is 
effectively "closer" to the resonance as the original pa- 
rameters would suggest. The current would then, in sec- 
ond order, be roughly given by the derivative of the first 
order term times the renormalization. 

The behaviour of the system at resonance (and its 
crossover to the Coulomb blockade regime) was also de- 
scribed in Ref. PJ7J] within the resonant tunneling ap- 
proximation for the two charge state model. Therefore, 
the expansion of the resonant tunneling formula up to 
second order yields Eqs. (|ll]) - (|l|) if we omit all terms 
with Ai and A_i. The integrals, then, become diver- 
gent and a cut-off (of the order of the charging energy) 
has to be introduced. In this Letter, however, we took 
into account all processes, and, therefore, no cut-off is 
needed. 

In Fig. g we show the second order contribution to the 
linear differential conductance G = dl/dV. Here and in 
the following we choose a symmetric coupling a\ = 0$. 
In Figs. I) and [| a comparison of the first order, the sum 
of the first and second order, and the resonant tunneling 
approximation (where the cut-off is adjusted as U = Eq) 
is displayed for the linear and nonlinear regime. The de- 
viation from the first order result (sequential tunneling) 
is significant and of the order 20%. The agreement with 
the resonant tunneling approximation provides a clear 
criterium for the choice of the bandwidth cut-off. Fur- 
thermore, and most importantly, it shows the existence 
of a parameter regime where renormalizations of Eq, do, 
and Ao by higher order charge states can be neglected al- 
though the current deviates significantly from the classi- 
cal result. We have checked the significance of third order 
terms ~ «g by using the resonant tunneling formula f@fil§ 
and exact results for the average charge in third order at 
zero temperature |15|| . For the parameter sets used in the 
figures, the deviations to the sum of first and second or- 
der terms were smaller than about 2%. Therefore, at not 
too low temperatures, second order perturbation theory 
is a good approximation even if the tunneling resistance 
approaches the quantum resistance. 

In Fig. p^ we present a comparison of our results with 
recent experiments [H . The temperature dependence of 
the Coulomb oscillations were measured for two different 
samples with conductances «o = 0.015 and ao — 0.063. 
For ao = 0.015, second order perturbation theory is suf- 
ficient and the results agree perfectly in the whole tem- 
perature and gate voltage regime. For ao — 0.063, third 
order terms start to become important, but the agree- 
ment is still reasonable. We emphasize, that only the 
bare values for ao and Eq have been used as they were 
determined unambiguously in the experiment. Further- 
more, the usage of the full resonant tunneling formula 



with the bare value of the charging energy would lead 
to a clear deviation from the experiment by about 10%. 
Thus, the inclusion of higher order charge states within 
second order perturbation theory, as presented in this 
Letter, is an important improvement of the theory. For 
«o = 0.063, the experiment is at the boarder line where 
third order terms and renormalization effects start to be- 
come important. We have checked this by using the reso- 
nant tunneling approximation with a renormalizcd value 
of the charging energy as estimated in the experiment. 
Again, the agreement is reasonable which indicates that 
higher order charge states lead to a renormalization of 
the charging energy. 

In conclusion we have presented a consistent calcula- 
tion of the tunneling current of the single electron tran- 
sistor up to second order perturbation theory. The ap- 
proach is free of any divergences and provides cut-off in- 
dependent results. At resonance we find new terms which 
are significant for experimentally realistic parameters. 
We have found a regime where second order perturbation 
theory can be used without any further renormalization 
of system parameters. A comparison with experiment 
shows good agreement. 
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FIG. 1. A diagram showing contributions to sequential tun- 
neling (S 'r and Sj ' ) and cotunneling (S ' and S '^~). 
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FIG. 2. The second-order contribution of the differential 

conductance G (2) = Y^Li ^ for T / E ° = °- 05 > a ° = °- 04 
and V = 0. 
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FIG. 3. The differential conductance for T/E c = 0.05, 
ao = 0.04 and V = 0: sequential tunneling, sequential plus 
cotunneling, and resonant tunneling approximation. 




FIG. 4. The differential conductance for T/E c = 0.01, 
ao = 0.02 and V/Ec = 0.2: sequential tunneling, sequen- 
tial plus cotunneling, and resonant tunneling approximation. 
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FIG. 5. Maximal and minimal linear conductance for 
E c = 1A7K and a = 0.015, and E c = IK and a = 0.063. 
The dots are experimental data from Ref. |l4j. 



